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We have simulated the possible O vacua states in heavy ion collisions. In a quench like scenario, 
random phases of the chiral fields were evolved in a zero temperature potential incorporating the 
breaking of Ua(X) symmetry. Initial random phases very quickly settles into oscillation around the 
values dictated by the potential. The simulation indicate that vacua states that can be populated 
in heavy ion collisions is a coherent superposition of a number of modes. 
QCD in the chiral limit (m u = m e i = m s = 0) possesses a Ua(X) symmetry. Spontaneous breaking of the symmetry 
require a neutral pseudoscalar Goldstone boson with mass less than v3iWj hi addition to pion itself. However no 
such Goldstone boson is seen in nature. The problem was resolved by the discovery of non-perturbative effcets that 
violates the extra J7a(1) symmetry, 't Hooft showed that because of the instanton solution of the Yang-Mills theory, 
the U(l) symmetry is not really a symmetry of the vacuum p|. Term containing the so called vacuum angle (0) 
which breaks P and CP symmetry can be added to the Lagrangian. QCD requires very small 9 ~ 10~ 9 , which 
explain the apparent P and CP symmetry in strong interaction. Dynamical breaking of the C/a(1) symmetry is also 
obtained in large N (color) limit of the SU(N) gauge theory In this approach the dominant fluctuations are not 

semi-classical but of quantum nature. 

Recently Kharzeev, Pisarski and Tytgat |5|,^| argued that in heavy ion collisions non trivial vacua states may 
be created. In the limit of large number of colors, the axial UaO-) symmetry of massless quarks may be restored at 
the the deconfining phase transition. As the system rolls back to confining phase, it may settles into a metastable 
state with non trivial 0. The idea is similar to the formation of Disoriented chiral condensate (DCC) j7j. In DCC 
space-time region is created where the chiral condensate points in a direction different from that of the ground state. 
Similarly, in vacua states space-time region with non-trivial may be created. 

After the suggestion of Kharzeev et al ||[| several authors have looked into various aspects of non-trivial vacua 
state that may be formed in heavy ion collisions j^-[nj]. Buckley et al 10, n umerically simulated the formation of 



vacua states, using the effective Lagrangian of Halperin and Zhitnitsky |12| . They assumed quench like scenario, 
rapid expansion of the fireball leave behind an effectively zero temperature region in the interior which is isolated 
from the true vacuum. Starting from an initial non-equilibrium state, they studied the evolution of phases of the 
chiral field. They saw formation of non-zero vacuum within a time scale of 10~ 23 sec. In their formulation they 
used a dissipative term with friction constant 7. 7=200 MeV was chosen, which may be rather large as the fields 
evolve essentially in a zero temperature potential, where dissipation will be small. They also choose to ignore the 
Fluctuation-dissipation theorem which require that dissipation be associated with fluctuations (noise). They did not 
include such a term. 

In the present paper we follow essentially the approach of Buckley et al |n],[ll| to investigate the possible creation 
of vacua states in heavy ion collision, with some important differences. We use the effective Lagrangian developed 
by Witten [[| to discuss the Ua(1) anamoly, and choose to omit the arbitrary dissipative term. In a quench like 
scenario, chiral phases are evolving in a zero temperature potential, where dissipative effects are supposed to be small. 
Also, neglecting dissipative term we are avoiding the problem of using a fluctuation term, non-inclusion of which will 
violate the fluctuation-dissipation theorem. We have also choose to use the proper time (r) and rapidity (Y) as the 
most appropriate coordinate system for heavy ion collisions. 

Effective non-linear sigma model which incorporates the breaking of symmetry can be written as Q, 

C = fUltr(d^d^U)+tr(M(U + U*)) - a(tr In U - 0) 2 ) (1) 

where U is a 3 x 3 unitary matrix with expansion, U = U${\ + i , }2 l t a TT a / fn + O^ 2 )), Uq being the vacuum expectation 
value of U, t a are the generators of U(3) (Trt a t b = S ab ) and 7r° are the nonet of Goldstone boson fields. M is the quark 
mass matrix, which is positive, real and diagonal. We denote the diagonal entries as /z 2 . They are the Goldstone 
boson sqared masses, if the anamoly term a (a/N in ref. j3|) were absent. Because M is diagonal U can be assumed 
to be diagonal , = e^'Sij. 
In terms of <^i's, the potential is, 

vfa) = ju- E cos & + °/ 2 (E & - Q ) 2 ) ( 2 ) 

It may be noted that as </>j arose from tr In U, it is defined modulo 2tt. In the present work we use fi 2 = 
(lUMeV) 2 , l4 = (161MeV0 2 ,^ 2 = (687) 2 and a = (A92MeV) 2 §,§. With these parameters, the mass matrix in 
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(0) can be diagonalised to obtain m° ~ 139 MeV, m v ~ 501 MeV and m v > ~ 983 MeV, close to their experimental 
values. 

Vacuum expectation values of the angles ^'s can be obtained from the minimisation condition, 

Mi sin ^ = o(0 - ^2<t>j) (3) 

Solutions of this non- linear coupled equations has been discussed in detail by Witten || . If the equation has only 
one solution then physics will be analytic as a function of 0. The solution vary periodically in with periodicity 
2ir ||]. Also this solution must be CP conserving whenever CP is a symmetry of the equation ||. However, it may 
happen that eq^ has more than one solution. Then the solution are not CP conserving, rather a CP transformation 
exchanges them. 

Existence of metastable states (0 vacua states) can be argued as follows: the vacuum expectation values depend 
on the ratio a//sf, which in the mean field theory decreases with temperature a//x| ~ (Td — T) 3 / 2 , Td being the 
decoupling temperature ||. Then as the system rolls back towards the chiral symmetry breaking, aj may be small 
enough to support metastable states. In ref. || it was shown that when a //if < .2467 there is a metastable solution, 
which is unstable in ir° direction unless a < a cri a cr /mu\ ~ .2403. 

Appropriate coordinates for heavy ion scattering are the proper time (r) and rapidity (Y) . Assuming boost invari- 
ance, equation of motion for the phases fa'B can be written as, 

ldfa d 2 <j>i d 2 (f>i v-^ 2 . ^ 

It is interesting to note that in this coordinate system, a dissipative term which decreases with (proper) time is 
effective. We have solved the coupled partial differential equations assuming a quench like scenario after the symmetry 
restoring phase transition. The rapid expansion of the high energy shell leaves behind an effectively zero temperature 
region, which is isolated from the true vacuum. The chirally symmetric fields then essentially evolves in the zero 
temperature potential. 

The set of partial equations (Q) were solved on a 64 2 lattice, with lattice spacing of a = lfm and time interval of 
a/10 fm. We also use periodic boundary condition. Solving eq.^ require field configuration at the initial time. We 
assume the initial time as Tj=l fm. Initial phases (<fii and <pi) were chosen according to following prescription, 

<fn = R N /(l + exp({r - r cr )/T) (5) 
<j> i =R N /{l + exp[{r-r cr )/T) (6) 

where Rn is a random number within the interval [— 27r/16, 27t/16]. We use r cr =20 fm and T—.5 fm. 

The phases were evolved according to the eq.^| for arbitrarily long time 50 fm in the zero temperature potential. 
The idea is to see whether nontrivial vacua emerge during the evolution. We have considered two situation, (i) 
0=0 and (ii)0=47r/16 within the radius r cr and zero beyond that radius. For 0=0 emergence of non-zero (j)'s will be 
signature of vacua. 

In fig. 1, the evolution of the spatially averaged phases are shown, 

/ <f>i(x,y)dxdy . 

fdxdy > z = u >"> s ( 7 ) 

For =0, the potential is minimised for <j> u — 4>d — <Ps — 0. The averaged fields also fluctuate around zero througout 
the evolution. This shows that chirally symmetric fields donot evolve into a non-trivial vacua. The situation is 
different for finite 0. For = 47r/16 the potential is minimised for </> n =0.502, 0^=0.243 and S =O.O13. Chirally 
symmetric fields very quickly reaches these values and oscillate around it. The oscillation continues for long, with 
little damping. As will be shown later, the continues oscillation of the spatially averaged phases indicate that a large 
number of modes contribute to the vacua state. 

We define a correlation function, 

CM = ( 8 ) 

such that the distance between the lattice points i and j is r. C (r) specifies how the three component (j) — (cf> u , phid, 4>s) 
at two lattice points are correlated. In fig. 2, we have shown the evolution of the correlation function for 0=0 and 
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0=0.785 at different times, r=l,6,ll,16 and 21 feu. Initially at t—1 fm, for both the cases, correlation length is about 
1 fm, the lattice spacing indicating the fact that initial fields were random, with no correlation. The correlation length 
increases at later times. But for 0=0, the increase is marginal. However for finite O, correlation length increases very 
rapidly, and assumes quite large values. Physically, for finite O all the fas try to align themselves in some direction 
(i.e. in the O direction) thereby giving very large correlation length even when they are separated by a large distance. 
This figure clearly demonstrated the possibility of parity odd bubbles formation. Initially random phases evolve such 
that they points in the same direction, forming a large parity odd bubble. 

To compare our results with that of Buckley et al |1(J , we have also studied the evolution of the different modes 
of the phases fa. At each alternate time step, we apply a fast fourier transform to the spatial data. The fourier 
transformed data are then integrated over the angles to obtain momentum distribution, 

1 f 2 " 

&{k) = — fa(k,6)d6,i = u,d,s (9) 
Jo 

Each mode was averaged over some narrow bin. In fig. 3, we have shown the evolution of 4> u for 0=0 (upper panel) 
and 0=0.785 (lower panel). For both the cases, the modes oscillate throughout the evolution. Initially amplitude of 
the oscillation decreases, but at later time it remains more or less same. This is indicative of the fact that dissipative 
term decreases with time. We also find that higher modes are suppressed compared to zero mode, but the suppression 
is not as large as obtained by Buckley et al j[(J[ll). Even after a large interval of time, higher modes donot become 
negligibly small. For 0=0, throughout the evolution, the zero mode as well as higher modes, oscillate about zero. 
But for finite value of O = 47r/16 the modes oscillate about some definite nonzero value. For (j)^ and fa, we have 
obtained qualitatively similar results. Indeed the evolution of different modes are in accordance with fig.l, where we 
have shown the evolution of spatial averaged phases. Spatially averaged phases were found to oscillate about some 
fixed value dictated by the minimum of the potential. Present simulation results are qualitatively different from the 
simulation results of Buckley et al |Tq]Tl|] . They found that for finite O, the modes quickly reaches the value dictated 
by the O vacuum. The difference is essentially due to the strong dissipative term in the equation of motion of Buckley 
et al JlO[ , which we have omitted. We have checked that if we include an arbitrary dissipative term, zero mode as 
well as higher modes quickly evolve into a constant value, dictated by the minimum of the potential. Also, in the 
simulation of Buckley et al [ |lO|]TT| , the O vacua is essentially consists of zero mode, while we find that in absence of 
dissipation, higher modes also contribute substantially. The O vacua is thus a coherent a superposition of a number 
of modes. 

Present simulation indicate that in heavy ion collision non-trivial ©-vacua state that is a coherent superposition of 
a number modes can be formed. What will be the signature of such a state. As such detecting O states are difficult. 
Kharzeev and Pisarski || estimated the P-odd observables are on the order of 10 -3 , a small effect. Also as discussed 
by Voloshin the so called signal of O states may be faked by "conventional" effect such as anisotropic flow etc. 
Whatever be the signal of the vacua states, with a large number of modes contributing, signal will be broadened, 
effectively diluting the detection probability. 

To summarise, we have simulated the non-linear sigma model incorporating the dynamical Ua{^) breaking. As- 
suming boost invariance equation of motion of the chiral phases were solved on a 64 2 lattice with lattice spacing a=l 
fm. It was shown that phases of chiral fields with finite O oscillate around the values dictated by the minimum of the 
potential. Correlation studies shows that initially uncorrelated phases very quickly develops large correlation length, 
indicating formation of a 'Parity Odd' bubbles. Fourier analysis of the modes indicate that the O vacua states are 
coherent superposition of a number of modes, which continue to oscillate with time. 
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FIG. 1. Evolution of spatially averaged 4>i, i=u,d and s for O=0 (solid line) and 0=0.785 (dashed line). 
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FIG. 2. Correlation function at different times for 0=0 (upper panel) and 0=0.785 (lower panel). 




FIG. 3. Evolution of the Fourier transformed <f> u (k) with time for different modes, k=6.2, 18.5, 30.9, 43.3 and 55.5 MeV. 
The solid line is for the lowest mode. 
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